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Diagnosis of Disease bv Det ermination of Electrical Network Properties 

of a Body Part 

Field of the invention 

This invention relates to a method for detecting and diagnosing disease 
states in body parts of living organisms by using a plurality of electrical 
impedance measurements. 

Bacl<around of the invention 

Measurement devices of various l<inds have been used to diagnose 
disease. For example, x-ray machines measure tissue density, ultrasound 
machines measure acoustic density, and thermal sensors measure 
differences in tissue heat generation and conduction; all such measurements 
can have diagnostic value. In addition to these devices, there exist devices 
that can measure electrical data, such as voltage between two points when a 
unit current is injected between other two points. If each point of the current 
injection Is in immediate proximity of one of the voltage measurement points 
the measured value is equal to the impedance of the body part between the 
cunrent injecting electrodes. 

Values of electrical impedance of various types of body tissue are well 
l^nown through in vivo studies on humans or from excised tissue made 
available following therapeutic surgical procedures. It is well documented that 
a decrease in electrical impedance occurs in the tissue as it undergoes 
cancerous changes. This finding is consistent over many animal species and 
tissue types, including, for example, human breast cancers. 

One technique for screening and diagnosing diseased states within the 
body using electrical Impedance is disclosed in U.S. Pat. No. 6,122,544, 
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which is incorporated herein by reference. In tliis patent, data are obtained 
from two anatomically homologous body regions, one of which may be 
affected by disease. One subset of the data so obtained is processed and 
analyzed by structuring the data values as elements of an impedance matrix. 
Such matrices can be further characterized by their eigenvalues and 
eigenvectors. These matrices and/or their eigenvalues and eigenvectors can 
be subjected to a pattern recognition process to match for known normal or 
abnormal (resulting from disease) matrix or eigenvalue and eigenvector 
patterns. The matrices and/or their eigenvalues and eigenvectors derived 
from each homologous body region can also be compared, respectively, to 
each other using various analytical methods and then subjected to criteria 
established for differentiating normal from diseased states. 

The published international patent application, PCT/CA01/01788, 
discloses a breast electrode array for diagnosing the presence of a disease 
state in a living organism, wherein the electrode array comprises a flexible 
body, a plurality of flexible arms extending from the body, and a plurality of 
electrodes provided by the flexible arms, where the electrodes are arranged 
on the arms to obtain impedance measurements between respective 
electrodes. In one embodiment, the plurality of flexible arms are spaced 
around the flexible body and provided with an electrode pair. In operation, the 
electrodes are selected so that the impedance data obtained will include 
elements of an nxn impedance matrix, plus other impedance values that are 
typically obtained with tetrapolar impedance measurements. In a preferred 
embodiment the differences between corresponding homologous impedance 
measurements in the two body parts are compared in a variety of ways that 
allow the calculation of metrics that can serve either as an indicator of the 
presence of disease or localize the disease to a specific breast quaidrant or 
sector. The impedance differences are also displayed graphically, for example 
in a frontal plane representation of the breast by partitioning the impedance 
differences into pixel elements throughout the plane. 
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However convenient this metliod is, obtaining impedance differences in 
this manner is simplistic. When current is injected, the flow of electricity 
through the body part is expected to take one or more complex, winding paths 
from the starting electrode to the final electrode. The conventional approach 
does not account for the internal current pathways that the electricity follows. 

Summarv of the invention 

The purpose of the procedure to be described Is to obtain a 
representation of a part of the human body in the form of an electric network 
or some equivalent to it. The usefulness of such a representation is due to the 
fact that physiological properties, such as alterations of structure due to a 
tumor, are generally associated with changes of electric conductivity. 
Therefore the resultant representation has diagnostic value for detection of 
anomalies. 

in particular, a system for diagnosing the possibility of disease in a 
body part is described. . The system includes a data acquisition module for 
acquiring the impedance matrix of the body part, as disclosed In U.S. Pat. No. 
6,122,544. The system further includes a network module for representing 
the body part by an electric networic of branches with initially unknown branch 
impedances, the network having external nodes corresponding to the location 
of the electrodes from the data acquisition module and internal nodes that 
cannot be accessed from the outside of the body part. The internal and 
external nodes are connected by current pathways. The system also includes 
an electrical properties module for determining electrical properties of the 
pathways using the measured electrical data, and a diagnosis module for 
utilizing the electrical properties to diagnose the possibility of disease in the 
body part. 
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The electrical properties module uses measured data and numerical 
techniques to determine the admittance (impedance) of each of the current 
pathways, thereby obtaining an admittance matrix for the particular network 
representation of the body part. The electrical properties module can likewise 
repeat these steps to obtain the admittance matrix associated with a 
homologous body part. For instance, the admittance matrices associated with 
the left and the right breast can be obtained and then compared by the 
diagnosis module to diagnose disease. 

Also described herein is a method for diagnosing the possibility of 
disease in a body part. The method includes measuring electrical data of the 
body part with a set of electrodes, and representing the body part by an 
electric network. The network has extemal nodes corresponding to the 
location of the electrodes and internal nodes, the internal and external nodes 
being connected by current pathways. The method also includes determining 
electrical properties of the pathways using the measured electrical data, and 
utilizing the electrical properties to diagnose the possibility of disease in the 
body part. 

Brief description of the drawings 

Figure 1 is a data flow diagram of the method for detecting and 
diagnosing the possibility of disease in a body part; 

Figure 2 is a sample network of nodes and current pathways, 
according to one embodiment of the present invention; 

Figure 3 is a data flow diagram of the electrical properties module of 
Figure 1 ; 
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Figure 4 is a data flow diagram of the diagnosis module of the 
diagnostic system of Figure 1; 

Figure 5 is a flow chart of the algorithm illustrating the method 
performed by the diagnostic system of Figure 1 to diagnose disease; and 

Figure 6 is a Current Pathway Difference (CRD) Plot for an actual 
subject in one embodiment of the present invention. 

Detailed description of the invention 

Figure 1 shows the outline of the proposed method for detecting and 
diagnosing disease in a body part, such as cancer in a breast. The method 
uses impedance measurements taken from the multi-channel impedance- 
measuring instrument 12 with the pair of electrode arrays 14 similar to the one 
described in PCT/CA0 1/0 1788, a network module 16, an electric properties 
module 18 and a diagnostic module 20. In one embodiment, the electrode 
array 14 includes Nc current Injection electrodes, and voltage measurement 
electrodes that are applied on the body part, each of the current injection 
electrodes being associated with the adjacent voltage measurement 
electrode. Two sets of measurements are performed, in the first, the 
impedance is measured between two voltage electrodes when the current is 
injected between associated current electrodes. The total number of 
independent current injections and related impedances in the first set is 
iVci =iVe (iVe-l)/2. In the second set of measurements the complete impedance 
matrix with respect to a base electrode Nb is obtained, where Nb can be taken 
to be Ne for example. For each of -1 current injections, i^, applied 
between current electrodes k and Nb , where k= l,2,...,Ne, k^Ns. voltages 
Vki between each of the voltage electrodes j = 1, 2,...,iVe and the base Nb are 
measured. In the second set of measurements there are Nc-1 current 
injections and for each current injection there are iv; - 1 voltage 
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measurements. The voltages thus obtained Tig are divided by the injected 
current ik. The resulting impedances can be organized in a so-called 
Impedance Matrix (IM) having elements Ztj = VylIj , as described in U.S. 
Patent No. 6,122,544. 



In the first set of measurements the impedance is measured wa times 
resulting in »ci impedance values, {Z*" , Zf Z;^ }, where Zf is the 

impedance measured between the voltage electrodes associated with the/* 
cunpent injection electrode pair when current is injected between associated 
current electrodes, as required in tetrapolar impedance measurement. 



From the second set of measurements the impedance matrix is 
obtained. In one embodiment of the invention where NB = Nothe following 
impedance matrix is obtained: 



"11 
^21 



^12 
^22 



^2(Ar,-l) 



(1) 



The choice of a base electrode is arbitrary. iVlultipie instances of iiVI for 
different base electrodes could be used to reduce the effect of measurement 
errors. 



Network module 16 includes hardware and/or software for representing 
the body part by a network. A number of different network topologies are 
possible with a given array 14. The topology to be used is an Input 
parameter. The network has extemal nodes corresponding to the location of 
electrodes in the electrode anray 14 and internal nodes. The extemal nodes lie 
on the perimeter of the network, while the internal nodes lie inside. A current 
pathway is a line segment that connects any two, and only two. nodes. 
Current pathways Intersect only at extemal nodes or internal nodes, and are 
the conduits through which current flows. The term "branch" will be used 
interchangeably with the term "current pathway." Each branch is associated 
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with two nodes and a branch admittance/impedance. A description of the 
network appears below in connection with Figure 2. 

The electrical properties module 18 includes software and/or hardware 
for determining electrical properties of the current pathways using the 
measured electrical data. For example, the electrical properties module 18 
can use the electrical data obtained by the multichannel impedance 
measuring instrument 12 to find the Impedance/admittance of each of the 
current pathways in the network as specified by the grid module 16. In the 
preferred embodiment of the invention, the network is represented by a 
conductance matrix. 

Diagnosis module 20 utilizes the determined electrical properties of the 
pathways to diagnose the possibility of disease in the body part. For 
example, diagnosis module 20 can compare the conductance matrix of the 
body part to an average matrix obtained from a population group, or to a 
conductance matrix obtained from a homologous body part, as described in 
more detail below. 

Figure 2 shows one possible network topology 50. The network 60 
includes eight external nodes 52, four internal nodes 54 and twenty-four 
current pathways 56. External nodes coincide with the location of the 
electrodes. All other nodes are called internal nodes. Internal nodes are not 
directly accessible and observable. The total number of external nodes is the 
same as the number of the current electrodes Ne, while the number of internal 
nodes depends on a chosen network topology. 

The electrical properties module 18 calculates electrical properties, 
such as the admittances of the pathways 56. The sample network 50 
represents the body part and the pathways represent equivalent paths where 
the current travels when injected into the body part. By using such a network, 
a more complete description can be obtained of the conductance properties of 
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the body part as compared to a network in which no internal nodes are 
included. In the example of Figure 2, the current pathways are line segments 
that intersect only at external nodes or internal nodes. 

Figure 3 shows the electrical properties module 18 of Figure 1. The 
electrical properties module 18 includes an admittance module 24, Zsame 
module 26, averaging module 28 and a conductance calculator 30. The 
admittance module 24 calculates an N^xN^ admittance matrix Y^^ from the 
impedance matrix Z^^^ provided by the multichannel Impedance measuring 
instrument 12 for the base electrode JV^, as known to those of ordinary skill. 

Zsame module 26 calculates the admittance matrix r^^ome from 

tetrapolar impedance measurements {Z^ . } as follows. Zsame is 

an upper triangular matrix of impedances whose non-zero elements can be 
represented as a vector, s. As mentioned above, the impedance matrix Z is 
the matrix obtained using one node as a base node, N^, which may be taken 
to be the last node. N. Zsame can be obtained fi-orh the elements of Z, z^j, 
according to 

Zsame = , - 2Z^^ + Zjj 
The upper triangular part of the impedance matrix can also be arranged in 
vector form, similar to the way s was obtained. The relationship between z 
and s is based on the last equation with indices set in the appropriate fashion 
for the vector form and is described using the matrix C that satisfies the 
equation: 

s = C*z 

The condition number of C is such that C"* can be derived with virtually no 
error. In this case, elements of the impedance matrix can be obtained from 
z^cr^s. The impedance matrix Z may now be obtained from z, making use of 
the fact that the impedance matrix is symmetric to obtain the lower triangular 
part. Once the impedance matrix in respect to N-th external node is obtained. 
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the admittance matrix can be easily obtained by inverting it and extending it 
by one row and column. 

Thie averaging module 28 calculates the average admittance matrix F 
from admittance matrices and Yj^^. In other embodiments, the average 
can be obtained by averaging Y^^^, Y^^ and other impedance matrices 
obtained with varying base electrodes. Another possibility Is to not average at 
all, and instead to use just one of the aforementioned impedance matrices. 

The conductance calculator 30 calculates the complete 
iN^ +Ni)x(N^ conductance matrix, G, for the body part by using the 

average admittance matrix F using Equation (4), as described below. The 
conductance matrix characterizes the conductance properties of the body 
part. Because the presence of some diseases, such as cancer, is known to 
alter the conductance of a body part, the conductance matrix possesses 
considerable diagnostic value. 

Figure 4 shows the diagnostic module 20 of Figure 1. The diagnostic 
module 20 includes hardware and/or software for comparing the branch 
conductances of the body part to conductances of corresponding branches of 
a typical body part obtained from a population group, making use of the fact 
that off-diagonal elements of the conductance matrix are equal to the 
respective branch conductances. A difference between corresponding branch 
conductances could indicate the presence of disease. In one embodiment of 
the invention, the diagnostic module 20 compares the conductance matrix G 
to a typical conductance matrix Grypicai, such as an average conductance 
matrix of an appropriate population group. Alternatively, or In addition, the 
diagnostic module 20 can compare the conductance matrix for the body part, 
such as the left breast Gl^ , to a conductance matrix for a homologous body 
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part, the right breast GM^t. Again, a difference between the two matrices 
could indicate the presence of disease. 

The diagnostic module 20 can calculate the absolute difference 
C^^Lea-G^RightJ or the relative difference [G^^-G^)^ ^ (Gua + GRigh,)] that is 

indicative of the possibility of disease in the body part or the homologous body 
part. Where there is a significant difference, further analysis can be 
performed to discern which of the homologous pair may be cancerous. For 
example, as described above, it is known that the electrical properties of 
cancerous tissue deviate from the norm in a predictable way. Thus, the body 
part having electrical properties more like those of a cancerous body part can 
be suspect. 

Figure 5 shows a flowchart that illustrates the main steps 70 utilized by 
system 10 to diagnose the possibility of disease in a body part. The first part 
of the procedure is preparatory and involves the acquisition of data for typical, 
normal body parts as follows: At step (71), the baseline body part is 
represented with a grid of current pathways. The grid can be two- 
dimensional, or three-dimensional. Next, at step (72), a number of healthy 
subjects are analyzed yielding a database (74) of impedances and associated 
branch admittances for a typical body part. We may repeat these steps to 
collect several typical sets of data depending on the size, body fat, or some 
other characteristic of the subject or the body part. This concludes the 
preparation part. The subject-specific part of the procedure follows: At step 
(76) a plurality of electrodes is applied to the body part, such as a breast and, 
at step (77), the electrodes measure the impedance of the body part between 
electrode pairs. At step (78), the branch impedances of each current pathway 
are calculated using electrical properties module (18). Subsequently, at step 
(70), the diagnostic module (20) is utilized to diagnose the possibility of 
disease in the body part. 
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Referring to Figs. 6A and 6B. sample results in the form of two Current 
Pathway Difference (CPD) plots are shown illustrating the value of the system 
and the method of the present invention In detecting breast cancer. In Figs. 
6A and 6B, the right and left breasts, respectively, of a female subject are 
represented as network plots. To generate these plots, each breast was 
represented by a network of lines representing current pathways. The branch 
admittances, as calculated by the diagnostic module 20, are plotted for each 
of the branches on the left breast element where the conductance of the 
branch on the left side is higher than the conductance of the analog branch on 
the right side and on the right breast where the conductance of the branch on 
the right side is higher than the conductance of the analog branch on the left 
side, othenvise just dashed lines are shown. The difference is plotted as a line 
with thickness proportional to the difference (the higher the difference the 
thicker the line). To obtain the branch admittances, use has been made of the 
fact known to those of ordinary skill that the branch admittances are non-zero 
off-diagonal elements of the admittance matrix. 

The CPD plot of Figure 6A indicates the presence of cancer in the 
upper inner quadrant of the right breast. (Biopsy confirmed a ductal 
carcinoma in this quadrant.) Observation of Figures 6A and 6B shows the 
preponderance of absolute differences in the diseased region of the right 
breast (most of the branches with lines), whereas Figure 6B shows most of its 
branches as dotted lines. 

In one embodiment of the invention the conductance calculator 30 from 
Figure 3 obtains the conductance matrix, G, that satisfies: 



where, if is the number of external nodes and N, is the number of internal 
nodes, then G„ is an iV^xiV, matrix, G,, is an N^xN, matrix and is an 




(3) 



(2) 
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N,xN, matrix. Here, v, is the vector of internal node potentials, is vector of 
external node potentials, is the vector of current injections into external 
nodes (this does not include the current through branches) and i, is the vector 
of currents that are injected into internal nodes (these are zero because these 
nodes are not accessible). 

The conductance calculator 30 obtains the conductance matrix by 
solving the conductance equation 

G„-G^'Gl=G^ (4) 

where G„ is the matrix of the real part of the measured admittances 
G„ =Ke^). Equation (4) is obtained from Equation (2) by eliminating the 
vector of internal potentials, v,. using v,==-G^^Glv^, and then substituting 

The conductance calculator 30 can solve Equation (4) using several 
methods. In the preferred embodiment of the invention Newton's method for 
solving a system of nonlinear equations is used. In the general case, there are 
more equations in (4) than unknowns, thus this set constitutes a linear least 
squares problem. This makes it Gauss-Newton. Its solution Is a well known 
iterative procedure. 

The only problem may occur when we do not have a good initial guess 
for conductances. Newton's method is guarantied to converge only when we 
are close to the solution. But when we are not, as in the case of a poor initial 
guess, then what we do is that we bring the "solution" closer to where we are. 
This may be done by a procedure known as a continuation method. The use 
of the continuation method assures robustness for the described 
methodology. 
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The Gauss-Newton method and a continuation method used in a 
preferred embodiment of the invention are now described in detail. 

Solution of Conductance Equation: 
The general form of the conductance equation (Eq. 4) is 

Ag) = (G^ - G,,GZ'Gl) - G„ = 0 (5) 

\N\\eref(g) is a function of the branch admittances to be an-anged in a vector^: 



g = 



See 
Sei 
Ste 



See 



Sei 



G^N,) 
\G^{:,N,)\ 





' Gl{:.,\) ' 


> Ste ~ 






G:i:,NJ 




' G,(-J) - 


> Sii = 




U 


Gui:>N,l 



(6) 



(A) Newton's Method for a system of nonlinear equations 
Let g = [^,], / = U";Nx N ,bea column vector of the conductances 
of all possible pathways for a given N. Let be a vector of estimated 
elements of G for the given network. Vector g has zero for each that 
corresponds to certain G(j},q) where nodes p and q are not connected. The 

Jacobianis = Substituting g„ in (5), /(^o)5tO is calculated. 

og 

A new vector g^"^ is calculated from: 
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j(^(«))Ag(")=-/(g('')) . 

For go close enough to g* (solution of the problem (5)), the iterative method 

described by equations (7) converges to g* (quadratically for a full set of 

equations). If g^ is not a good initial guess for g*, system (7) may not 

converge. In that case, the continuation method described below is applied. 
If the number of equations is greater than the number of unknowns, a least 
squares problem is being solved using the Gauss-Newton method. In that 
case, the convergence is not quadratic. 

(B) Continuation Method 

For a function f :R"-^ R"' , a numerical method, sensitive to the initial 
guess, can be used to solve the equation /(g) = o , wheregG R" . If the guess 
is not sufficiently close to the solution g* , the numerical method might not 
converge- To obtain a sufficiently close initial vector g^^^ for problem 

/(g) = 0 , an augmented function fig(d\0) : R"^^ R"* can be considered 
that is defined as: 

hgie\e)^{\-eyf{gm^ef{gm , (8) 

and the following equivalent problem solved: 

7(g(0),0) = O, 0e[O,l] (9) 

where: 

f(g):R" -^R"* is a function of g, such that equation /'(g) = 0 is easy to 

solve, and 0 is a real number. 

The important properties of the new function are: 
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As), 
fis)> 



0€]O,1[ 



0=0 



0=1 



For 0 = 0, the equation (9) becomes 7(^) = 0. which is easy to solve by its 
construction. 

For 0 = 1, the problem (9) Is identical to the original problem /(^) = 0 . 
For 0e]O,l[ , the problem (9) becomes 



The following method takes the solution of the problem^ = o as the initial 
guess g^^^for the close problem (l-0)g(0)+0 /(g(0)) = O . where is 6 =0+86 , 
and 50«1 . The underlying assumption Is that the solution for the problem 

7(g(0),0) is close enough to ^(0 + 50) for the problem fig(0+Se%e + 6d) so 
that the numerical algorithm converges tog* (0+50) . 

To summarize, the continuation method includes a) solving the 

problem /(^(0),0) = O for 0 = 0, b) taking the solution g*(0) as the initial 
guess g-^"^ (0 + 50) for 

f(g(0 + Se\e + 50) = (1 - (0 + 50)) 7(g(0 + 50)) + (0 + 50) /(g(0 + 50)) . (11) 

where 50 is a small step, c) attempting to solve this problem using a Gauss- 
Newton algorithm, d) if the Gauss-Newton algorithm converges, setting 
0=0+50, Sd = a de, where a>l is an acceleration coefficient, and if it 
does not converge, setting 39=89/2, and e) repeating step b) until 0=1. 



(i-e)/(g(0))+0/(g(0))=o 



(10) 
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These methods can be used to solve the conductance equation: 

fis) = (.Gee - G,,G-^Gl) -G„ (12) 
The matrix of finite differences \s 6f = ]6f^ bf^ . . • ] , where 

From the above equation, 

4/* -^G^G;'g^ -G^AG;,'g^ -G^G^'Ag^ (13) 

Matrix AGj' can be calculated from the equation g;^g„ =/as follows 

G-'G„=I |A 

AG-'G„+G-'AG„=0 \G-' 



AG;»=-G,7'AG„G;' (14) 
Substituting (14) in (13) leads to: 

Af, =A5«*-AG,,(G-^g^) + (G,,G«')AG,(G-'5^,)-(G.,G-')A^,,, (16) 

The next step is to reorder the matrix products in (15). Each AG„ should 
come as a vector Ag«to the right hand side of the conresponding matrix 
product. Aff*. is a column vector of all columns of AG^ as indicated by (6). 
By applying the rule 
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AXB <=> AAx (16) 

Equation (16) becomes 

^fk = ^eek - M^^g^i + M^^gi, - M^Ag,,, (1 7) 

where =G^GJ^ , M^^ =G^G;U^ . and M,^ and are calculated by 
applying (16): 

AG^iG;:'g^)=M^Ag^ 
AG„(G;^g^) = A,Ag„ 



Equation (17) can be written as: 



(18) 



where = [/ -m^ ji/^j^ -aTj] is the Jacobian matrix corresponding to the 
column of /(g^"^) in equation (7). 



Solving the system 



is equivalent to solving equation (19) 



I 



I 



(19) 



To simplify equation (19), the vector of admittances may be written as 

gg = T,-g^ (20) 
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gg- 



18- 



8l, 



In (20) is a vector of the current pathways in the selected network, 
while is an incidence matrix defined by the network so that equation (20) 
holds. 



Substituting (20) in (19) results in the reduced system (21): 



(n) 



where 



A/;' 



0 I -M<f' I M^I' i -Mf 
I -<> I I 
1 i I : I 



(21) 



0 



and i^t^a, 's a vector of current pathways admittance corrections. The new 
admittance vector g^'^ is calculated from the equation 

Sbrandi S branch ^^Z,^ (22) 



Similar to (20), vector g^"*^^ can be calculated from g^^^ as 

„(»>+>) _ T . (fi.-*^) 

where 7* is a matrix defined by the chosen network so that the equation (23) 
holds. 
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There are other methods, besides the method just described, which 
can be used by the conductance calculator 30 for obtaining th^ branch 
conductances and subsequently the conductance matrix. 

Once the conductance matrix calculator 30 finds the conductance 
matrix for the body part, using any of the methods described above, the 
aforementioned steps can be repeated to obtain the conductance matrix of 
the homologous body part. The diagnosis module 20 can then compare the 
conductance matrix for the body part and the conductance matrix for the 
homologous body part by using several comparison methods. For example, 
the norm of the difference of these two matrices can be computed, and if it is 
greater than some threshold, then further analysis can be performed as this 
difference may signal the presence of disease. 

Several computer systems can be used to implement the method for 
diagnosing disease in a body part. The computer system can include a 
monitor for displaying parts or the whole conductance matrix, or for displaying 
the difference between the conductance matrix for the body part and the 
conductance matrix for the homologous body part using one of several visual 
methods. In one embodiment, the method can be Implemented on a 2 GHz 
Pentium™ 4 system with 512 MB RAM. 

It should be understood that various modifications could be made to 
the embodiments described and illustrated herein, without departing from the 
present invention, the scope of which Is defined in the appended claims. For 
example, although emphasis has been placed on describing a system for 
diagnosing breast cancer, the principles of the present invention can also be 
advantageously applied to other diseases of other body parts. 

Additionally, in the above analysis, emphasis was placed on the real 
conductance matrix G, which was obtained by solving the conductance 
equation (4); however, it should be understood that an analogous analysis in 
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the complex domain could also be carried out by obtaining instead the 
complex admittance matrix Y, given by 



Y = 



by solving the admittance equation Y^^-Y^^%=Y, where Y is the average 
and y^teme- The complex admittance matrix may then be used to 
diagnose disease in a manner analogous to the way its real part was used 
above. In the last equation, the average impedance matrix need not be used; 
Instead, one of the aforementioned impedance matrices, Yj^^ or Y^^ can be 
used as the right hand side of this last equation. In this embodiment, where Y 
Instead of G is computed, analogous steps are taken to those described 
above, mutatis mutandis (for example, replacing the conductance calculator 
30 by an impedance calculator 30). 



